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LONG-TERM  GOALS 

Modeling  and  predicting  ocean  optical  properties  requires  linking  optical  properties  with  the  physical, 
chemical,  and  biological  processes  in  the  upper  ocean.  Our  long-term  goal  is  to  incorporate  optical 
processes  into  coupled  physical-biological  models  for  both  open  ocean  and  coastal  waters,  develop  and 
improve  integrated  ocean  forecasting  systems,  including  prediction  of  ocean  optical  properties. 

OBJECTIVES 

1)  To  improve  performance  of  the  coupled  physical-biological  model,  which  is  based  on  the  Navy 
Coastal  Ocean  Model  (NCOM)  for  the  California  Current  System  and  Regional  Ocean  Model  System 
(ROMS)  for  the  Pacific  Ocean; 

2)  To  incorporate  optical  variables  into  the  improved  coupled  3D  physical-biological  model  for  the 
Pacific  Ocean  and  California  Current  System; 

3)  To  evaluate  physical-biological-optical  models  with  remote  sensing  and  available  in  situ 
observations; 

4)  To  use  these  variables  to  drive  a  radiative  transfer  model  (EcoLight)  that  will  simulate  and  predict 
the  subsurface  light  field  as  well  as  the  ocean’s  color. 

APPROACH 

To  achieve  the  first  objective,  we  have  conducted  a  series  of  3D  physial-biological  model  simulations 
for  the  Pacific  Ocean  and  California  Current  System.  We  have  incorporated  the  Carbon,  Silicate,  and 
Nitrogen  Ecosystem  (CoSiNE)  model  into  both  the  Navy  Coastal  Ocean  Model  (NCOM)  and  Regional 
Ocean  Model  System  (ROMS).  The  CoSiNE  model  (Chai  et  ah,  2002;  2003,  and  2007)  was  developed 
originally  for  the  equatorial  and  Pacific  Ocean.  During  the  past  several  years,  we  have  implemented 
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the  CoSiNE  model  and  improved  the  performance,  and  investigated  the  seasonal  and  interannual 
variation  of  biological  productivity  in  response  to  both  local  and  remote  climate  forcing  (Chai  et  ah, 
2009;  Liu  and  Chai,  2009a,  b). 


To  simulate  optical  variables  with  the  improved  coupled  3D  physical-biological  model,  we  incorporate 
spectrally-resolved  inherent  optical  properties  (IOPs)  into  the  existing  ROMS-CoSiNE  model,  which 
haven’t  included  the  optical  processes  and  in-water  radiative  transfer  model  yet  (see  Figure  1).  We 
have  successfully  incorporated  an  optical  module  along  with  a  commercially  available  radiative 
transfer  model  (EcoLight)  into  a  one-dimensional  ecosystem  model  which  is  flexible  enough  that  can 
be  applied  to  three-dimensional  circulation  model  (Fujii  et  ah,  2007).  We  are  collaborating  with  Dr. 
Curt  Mobley  at  Sequoia  Scientific  to  implement  the  updated  version  of  EcoLight  into  the  ROMS- 
CoSiNE  model.  To  start,  we  will  use  an  idealized  test  bed  framework  for  3D  ROMS-CoSiNE  model. 
Recently,  the  CoSiNE  model  has  been  fully  incorporated  into  the  standard  ROMS  code,  and  the 
CoSiNE  model  results  have  been  evaluated  and  compared  with  several  other  ecosystems  models  that 
already  imbedded  into  the  ROMS.  Later  phase  of  our  collaboration  with  Dr.  Mobley  will  imbed 
EcoLight  into  the  ROMS-CoSiNE  model  for  application  to  the  Pacific  Ocean  and  California  Current 
System. 


The  optical  model  explicitly  represents  spectrally-resolved  IOPs  including  absorption,  scattering,  and 
attenuation.  The  optical  model  requires  input  from  ROMS-CoSiNE  model,  such  as  phytoplankton 
associated  chlorophyll  and  non-algal  particles  (NAP),  and  it  provides  outputs  such  as  phytoplankton 
absorption  coefficient,  detritus  absorption  coefficient,  particulate  backscattering  coefficient.  These 
optical  products  are  all  spectrally-resolved  and  can  be  compared  with  satellite  and  in  situ  data  directly. 
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Figure  1:  Overall  approach  for  incorporating  optical  properties  into  the  coupled  three-dimensional 
ROMS-CoSiNE  model.  PI,  P2,  PS,  POC,  PIC,  and  CDOC  denote  pico-plankton,  diatom, 
coccolithophorids,  particulate  organic  carbon,  particulate  inorganic  carbon  and  color  dissolved 

organic  carbon,  respectively. 
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Since  the  last  annual  report,  we  started  developing  an  optical  model  that  is  fully  compatible  with 
ROMS-CoSiNE  and  NCOM-CoSiNE  model,  which  explicitly  represents  spectrally-resolved  IOPs. 
However,  the  weakness  of  the  early  model  is  that  the  color  dissolved  organic  matter  (CDOM)  was  set 
constant  at  that  time  due  to  its  complicated  optical  and  biogeochemical  properties.  CDOM  plays  an 
important  role  in  absorption  spectra  both  in  the  open  ocean  and  coastal  regions.  In  order  to  improve  the 
existing  model  and  try  to  simulate  CDOM  dynamics,  we  have  modified  CoSiNE  model  and 
incorporated  color  dissolved  organic  carbon  (CDOC)  as  a  state  variable,  which  is  expressed  as: 

BIOLOGY(CDOC)  =  BIOLOGY{P )  +  BIOLOGY(Z)  +  BIOLOGY (B)  -  UVlight  (1) 

Here,  CDOC  is  mainly  produced  by  phytoplankton  (P),  zooplankton  (Z)  and  bacteria  (B)  through  their 
biological  processes  such  as  grazing  loss,  excretion,  and  mortality  etc.  CDOC  is  destroyed  by  the  UV 
light  only,  which  is  a  function  of  the  light  attenuation  coefficient  (kd)  at  300  mn: 


UVIight(z)  =  aCD0C(4 10,z)  xi?x  exp[-|  k,(300,z)dz]  (2) 

0 

where  R  is  the  coefficient,  acDOc(410)  is  the  CDOC  absorption  coefficient  which  can  be  related  to 
CDOC  concentration  by 


acooc (4 10)  =  a*coc(410)  x  CDOC  (3) 

where  a*cDOc(410)  is  the  specific  CDOC  absorption  coefficient  at  410  mn.  Then  for  the  whole  visible 
spectra,  CDOC  is  modeled  as 


amoM)  =  <W(410)  x  exp {-0.0 145  x  (A -410)}  (4) 

In  order  to  better  simulate  water  optical  properties,  we  also  included  coccolithophorids  as  a  state 
variable  in  the  CoSiNE  model.  Phytoplankton  absorption  spectra  is  then  modeled  as: 

a^(A,z)  =  «’,(/ l,z)  x  chl\{z )  +  a*2(/l,z)  x  chl2(z )  +  a*3(/ t,z)  x  ch/3(z)  (5) 

where  chll,  chl2,  and  chl3  is  pico-plankton,  diatom,  coccolithophorids,  respectively,  a*®  is 
chlorophyll-specific  absorption  coefficient,  which  is  modeled  as  When  taking  into  account 
phytoplankton  photo-adaptive  stage: 
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where  3.  cpl (high  light)  (^)?  ^  tpl(low  light)  (^)?  &  (p2(high  light)  (^)?  3  (p2 (low  light)  (^)  ^I*e  chlOFOphyll-SpeciflC 

3bsorption  coefficient  3t  high  3nd  low  light  conditions  by  ench  phytoplnnkton,  respectively.  0cmin  and 
0  max  are  the  minimum  and  maximum  phytoplanktonic  chlorophyll  to  carbon  ratios.  9si  and  0S2  are 
chlorophyll  to  carbon  ratios  for  pico-plankton  and  diatoms,  respectively.  Note  that,  we  assume 
coccolithophorids  do  not  have  photo-adaptive  effect  in  the  current  model. 
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In  the  optical  model,  assuming  no  contribution  of  CDOM  to  backscattering,  backscattering  coefficient 
is  determined  as  the  sum  of  seawater,  algae  and  co-varying  particles.  Algae  and  co-varying  particles 
are  expressed  as  a  function  of  particulate  organic  carbon  concentration  (POC)  and  particulate  inorganic 
carbon  (PIC)  of  small  and  large  particles,  which  consist  of  pico-plankton,  diatoms,  coccolithophorids 
and  associated  NAP,  respectively.  Backscattering  by  small  and  large  POC  (POC  1  and  POC2)  and  PIC 
are  formulated  as  follows  (Stramski  et  al.,  1999;  Balch  et  al.,  1996): 


,  ,  „  ,  ,  POC\(z)  '  ,  A  x_05 

^  ’Z  _  476935.8  X  510 

(8) 

bbl(A,z)  =  (P°C2^)^59 
b”2  ’  v  17069.0 

(9) 

bbp3(A,z)  =  (1.6  x  PIC(z)  -0.0036)  x  (^)!35 

(10) 

The  total  backscattering  coefficient  is  calculated 

bb(A,z)  =  0.5  b  J  A)  +  bhp](X,z)  +  bbp2(A,z)  +  bhJA,z)  +  bblig  (11) 

where  bsw(A)  is  seawater  scattering  coefficient  with  the  correction  for  salts  (Boss  and  Pegau,  2001). 
bbbg  is  the  background  backscattering  coefficient  the  implicitly  reflects  contributions  by  non- 
phytoplankton-covarying  bacteria  and  other  particles  (Behrenfeld  et  al.,  2005). 

WORK  COMPLETED 

This  report  summarizes  modeling  activities  between  1  October  2008  and  30  September  2009. 

We  have  completed  the  following  modeling  tasks: 

1 .  Evaluated  the  3D  ROMS-CoSiNE  model  performance  for  the  Pacific  Ocean  and  compared  the 
model  results  with  available  in  situ  observations  on  nutrients,  phytoplankton,  and  carbon  fluxes  (Chai 
et  al.,  2009,  Liu  and  Chai,  2009a, b,  Bidegare  et  al.,  2009;  Palacz  et  al.,  in  press). 

2.  Linked  the  bio-optical  module,  based  on  ID  configuration  (Fujii  et  al.,  2007),  with  the  ROMS- 
CoSiNE  model,  and  incorporated  CDOC  into  the  CoSiNE  and  bio-optical  module  for  the  Pacific 
Ocean; 

3.  Evaluated  the  ROMS-CoSiNE-Optical  model  results  with  remote  sensing  derived  IOPs  for  the 
Pacific  Ocean; 

4.  Incorporated  bio-optical  module  and  CoSiNE  into  the  Navy  Coastal  Ocean  Model  (NCOM)  for 
the  California  Current  System  (CCS); 

5.  Started  collaboration  with  Dr.  Curt  Mobley  to  incorporate  EcoLight  into  an  idealized  ROMS- 
CoSiNE  model. 
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RESULTS 


We  have  conducted  a  series  of  the  ROMS-CoSiNE  model,  without  the  optical  component,  for  the 
Pacific  Ocean  for  the  period  of  1990  to  2008.  For  doing  so,  we  can  evaluate  the  ROMS-CoSiNE  model 
results  with  available  observations,  and  then  improve  the  CoSiNE  model  perfonnance.  For  example, 
we  compared  the  modeled  nutrients,  phytoplankton  biomass,  and  carbon  fluxes  with  the  South  East 
Asia  Time-Series  Station  (SEATS)  in  the  South  China  Sea  (Liu  and  Chai,  2009a,  Chai  et  ah,  2009). 

We  evaluated  the  ROMS-CoSiNE  model  results  with  the  Hawaii  Ocean  Time-Series  (HOTS)  and 
addressed  how  ecosystem  responded  to  Pacific  decadal  climate  changes  (Polovina  et  ah,  2008, 

Bidigare  et  ah,  2009).  We  discussed  how  the  ROMS-CoSiNE  model  responds  to  1997-98  El  Nino 
event  in  the  Japan/East  Sea  (Liu  and  Chai,  2009b).  We  calculated  the  nutrient  transports  and  their 
impact  on  phytoplankton  productivity  in  the  equatorial  Pacific  (Palacz  et  ah,  in  press).  Since  these 
ROMS-CoSiNE  model  results  have  been  published  in  the  peer-reviewed  journals,  we  are  not  including 
these  results  in  this  report.  Here,  we  report  the  results  related  to  the  optical  model  development  and 
model  evaluation  with  remote  sensing  observations. 

Performances  of  the  coupled  ROMS-CoSINE-Optics  model  are  evaluated  by  comparing  the  model 
results  with  synchronous  MODIS  AQUA  measured  data  provided  by  NRL.  There  are  a  lot  of 
algorithms  being  used  to  calculate  ocean  optical  properties  from  satellite  data.  What  we  used  for  the 
model-satellite  data  comparison  is  the  one  called  quasi-analytical  algorithm  (QAA),  which  is  a 
promising  algorithm  for  deriving  inherent  optical  properties  from  ocean  color.  Figure  2  shows  the  data 
comparison  between  model  and  satellite  for  January  2006.  Spatial  averaged  (domain  shown  in  Fig.  2) 
total  absorption  a  is  0.035  (m'1)  at  440  mn  from  the  model  results,  and  0.03  (m'1)  at  443  mn  for  the 
MODIS  data  (based  on  QAA);  CDOM  absorption  coefficient  is  0.025  (m'1)  at  410  mn  for  the  model 
and  0.018  (m'1)  at  412  mn  for  the  MODIS  data;  the  particulate  backscattering  coefficient  bbP  is  0.003 
(m'1)  at  550  mn  for  the  model  and  0.002  (m'1)  at  555  mn  for  the  MODIS  data.  Overall,  the  optical 
model  produced  the  IOPs  tend  to  agree  well  with  the  MODIS  derived  results  in  terms  of  the  magnitude 
and  spatial  patterns.  We  continue  to  adjust  the  parameters  values  and  model  structure  to  improve  the 
model  performance  based  on  the  remote  sensing  derived  IOPs. 
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Figure  2:  Monthly  (January  2006)  averaged  comparison  between  the  modeled  and  QAA  derived 
chlorophyll  concentration,  particulate  organic  carbon,  CD OM  absorption  coefficient  at  410  nm, 
total  absorption  coefficient  at  440  nm,  and  particulate  backscattering  coefficient  at  550  nm  bhP(550). 
The  left  column  is  the  ROMS-CoSINE-optics  model  outputs,  the  central  column  is  MODIS  QAA 
products,  and  the  right  column  is  the  statistical  comparison.  The  colors  in  the  statistical  analysis 

( the  right  column)  denote  the  density  distribution. 
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Beside  the  incorporation  of  the  CoSiNE  model  into  the  ROMS,  we  have  been  collaborating  with  the 
Dynamics  of  Coupled  Processes  Section  (led  by  Dr.  Igor  Shulman)  at  NRL  to  implement  our  CoSINE 
model  into  the  Navy  Coastal  Ocean  Model  (NCOM)  for  the  California  Current  System  (CCS).  The 
coupled  physical  and  biological  model  assimilates  remote  sensing  information,  and  produces  near-real 
time  simulations  on  regular  basis.  The  physical-biological  model  results  can  be  viewed  and  obtained 
directly  from  the  NRL  website  (http://www7320.nrlssc. navy.mil/biospace/ModelingApproach.html). 
Penta  et  al.  (2008)  used  the  NCOM-CoSiNE  model  results  to  link  with  a  new  model  for  simulating 
underwater  light  attenuation,  which  improves  the  model  performance  for  the  Monterey  Bay.  The 
simulation  of  a  deep  chlorophyll  maximum  demonstrates  the  effect  of  the  new  light  model  at  depth 
(Penta  et  al.,  2008).  We  continue  our  collaboration  with  NRL  group,  and  will  send  our  modified 
version  of  CoSiNE-Optical  codes  at  regular  bases. 

IMPACT/APPLICATIONS 

Incorporating  ocean  optical  processes  into  coupled  physical-biological  models  will  enable  us  to 
simulate  and  forecast  optical  properties  in  coastal  waters.  With  demonstration  of  some  initial  successes 
of  developing  physical-biological-optical  modeling  and  data  assimilation  capability  for  the  Pacific 
Ocean  and  California  Current  System,  we  should  be  able  to  start  the  development  of  an  end-to-end 
ocean  forecasting  system.  Such  modeling  system  would  be  a  powerful  tool  to  design  the  adaptive 
sampling  strategy  and  would  be  an  essential  component  of  future  field  experiments  both  in  and  outside 
Monterey  Bay. 

TRANSITIONS 

We  have  been  working  with  Drs.  Igor  Shulman,  Bred  Penta,  and  Sergio  Derada  at  the  Naval  Research 
Laboratory  (NRL).  The  ecosystem  model  code  has  been  transferred  to  the  Dynamics  of  Coupled 
Processes  group  (led  by  Dr.  Igor  Shulman)  at  the  NRL.  Once  the  optical  module  is  fully  tested  with  the 
CoSiNE,  then  the  optical  module  along  with  the  updated  version  of  CoSiNE  model  will  be  transferred 
to  NRL. 

RELATED  PROJECTS 

This  project  has  strong  collaboration  with  other  ONR  supported  projects.  Besides  working  closely  with 
the  modeling  group  at  the  NRL  and  their  BioSpace  project,  we  are  collaborating  with  Dr.  Curtis 
Mobley  of  Sequoia  Scientific  on  improving  the  link  between  the  radiative  transfer  model  (EcoLight) 
within  the  ROMS-CoSiNE-Optics.  Recently,  Dr.  Mobley  has  obtained  a  grant  from  ONR  to  improve 
the  efficiency  of  the  EcoLight  calculation  in  the  3D  framework.  We  are  also  collaborating  with 
scientists  at  the  Monterey  Bay  Aquarium  Research  Institute  (MBARI)  to  use  the  observational  data  for 
the  region.  Dr.  Yi  Chao  at  JPL  has  been  collaborating  with  us  about  implementing  the  CoSiNE  into  the 
ROMS  for  the  Pacific  Ocean  and  the  CCS. 
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